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Abstract 

All the next-to-leading order results on Altarelli-Parisi splitting functions have been 
obtained in the literature either by using the operator product expansion method or 
by making use of the Curci Furmanski Petronzio (CFP) formalism in conjunction 
with light-like axial gauge, principal value (PV) prescription and dimensional re- 
gularization. In this paper we present the calculation of some non-singlet two-loop 
anomalous dimensions within the CFP formalism using light-cone axial gauge with 
Mandelstam-Leibbrandt (ML) prescription. We make a detailed comparison between 
the intermediate results given by the (PV) versus the (ML) method. We point out 
that the (ML) method is completely consistent and avoids the "phenomenological 
rules" used in the case of (PV) regularization. 
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1 Introduction 



At present and future high energy colliders a number of important hard scattering cross 
sections can be measured with high accuracy. The satisfactory description of these pre- 
cision data requires the evaluation of next-to-leading order (NLO) and in few cases even 
next-to-next-to leading order (NNLO) corrections in perturbative QCD. A number of 
NLO corrections are by now available in the literature (see e.g. [[[J and references therein). 
NNLO corrections, however, could be evaluated only in few cases. The special tech- 
nique developed for the calculation of the NNNLO corrections to the total cross section 
of electron positron annihilation into hadrons has recently been applied by Larin and 
Vermaseren || to get the NNLO and NNNLO anomalous dimensions for those operators 
which contribute to the Bjorken sum rule and the Gross-Llewellyn Smith sum rule. 

In the case of hard processes with two initial hadrons, however, no complete NNLO 
result is available. Although the NNLO coefficient functions of the Drell-Yan process 
have been obtained in ref. M, the phenomenological application of this result requires 
the calculation of the NNLO Altarelli-Parisi splitting functions (three-loop anomalous 
dimensions) as well. It seems that significant technical development has to be achieved 
before this calculation can be carried out. 

Two rather different methods , [§, have been used for the evaluation of the spin in- 
dependent two-loop anomalous dimensions. Both methods recently also have been applied 
to the successful calculation of the NLO corrections of the spin dependent Altarelli-Parisi 
splitting functions || ||.(] In the case of the first method || (OPE) one has to evaluate 
the overall ultraviolet divergences of twist-two local operator insertions. The results are 
obtained in moment space and the calculation can be carried out in Feynman gauge. Un- 
fortunately, the number of the operator insertions increases very rapidly in higher orders 
and the treatment of operator mixing becomes increasingly cumbersome The second 
method || [7[] (CFP) is based on the observation that in axial gauge the two-particle ir- 
reducible kernel of the ladder diagrams is finite. Therefore, using renormalization group 
properties, the anomalous dimensions are given by some projection of this kernel in agree- 
ment with the factorization theorem of mass singularities fl3 |. The corrections to the 
Altarelli-Parisi splitting functions are obtained directly in configuration space by evaluat- 
ing few Feynman diagrams. A detailed documentation of the NLO calculation using the 

2 A third method proposed in ref. Jlpfl makes use of the the operator definition of the parton number 
densities and requires the evaluation of Feynam diagrams with attached eikonal factors. The method has 
not yet been tested in higher order calculations. 
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CFP method appeared recently in ref. [ I2[ . A theoretically questionable feature of the 
(CFP) method is the use of light-like axial gauge with principal value prescription. It is 
not clear whether the PV prescription remains valid in higher orders since standard quan- 
tization procedures in light-like axial gauge lead to the so-called Mandelstam-Leibbrandt 



(ML) prescription [[14], Another technically unpleasant feature of the method is that 
the individual diagrams contributing to the finite 2PI kernel have soft and collinear sin- 
gularities (which cancel in the sum). The actual algebraic complexity of the CFP method 
nevertheless is significantly smaller than the one of the OPE method, motivating the 
further study of its technical aspects. In particular, it is of interest to investigate the 
CFP method with ML prescription. A first attempt already appeared in the literature: 
The one-loop Altarelli-Parisi splitting functions have been calculated by Bassetto [TBI . In 
the NLO calculation with ML prescription, the treatment of the ultraviolet part becomes 
more consistent. This is a clear principal advantage. However, the ML prescription is 
not without difficulties. Firstly, using two light-like vectors the integrals become more 
complicated. Secondly, the singular 1/nq factor appearing in the expression for the axial 
gauge propagator is regularized with a +ie prescription which results in new so-called 
axial ghost contributions. 

In this paper we apply the CFP method with ML prescription to the calculation of 
non-singlet two-loop anomalous dimensions. In particular, we calculate all terms which 
are proportional to C 2 F . An interesting new feature of our calculation is the evaluation of 
the axial ghost contributions. In Section 2 we give a short review of the CFP formalism, 
in Section 3 the PV and ML prescriptions are discussed. Section 4 describes details of 
the calculation, treating separately the virtual and the real contributions. In Section 5 
we make a crtitical comparison between the results obtained by using the PV and ML 
prescription and draw our conclusions from this discussion. An Appendix contains virtual 
integrals and phase space parametrizations. 



2 Framework of the calculation 

According to the factorization theorem, simple hard scattering cross sections can be writ- 
ten in perturbative QCD as convolutions of finite hard scattering cross sections and sin- 
gular parton number densities 



a{%, a 5 (/i 2 ), e) = a{%, a s (^)) ® r(a 5 (/i 2 ), e) (1) 
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where the symbol <8> indicates the convolution over longitudinal momentum fractions which 
in moment space is reduced to a simple product 



cr A r(^-,a 5 ( y u 2 ),e) = 6-Ar(^-,a5(/i 2 ))rAr(a5(/i 2 ),e) (2) 



where 



f N = f 1 dxx N ~ l f(x). (3) 

JO 

Our discussion will be valid for the non-singlet parton densities, but it can be easily 
generalized to the singlet case (matrix problem). The physical cross section an must 
be independent from /z 2 , therefore the deviation of the physical dimension of a from its 
canonical value (anomalous dimension) 



d Q 2 

-f N (as(^ 2 )) = --7z ^lna N ( — ,a s (^ 2 )) (4) 

d In 



is related to the moments of parton number densities 

7Ar(as(/i 2 ) = — lnrjv(as(/i 2 ),e) = f3(a s , e)-^— In T N (a s (/J, 2 ), e) (5) 
dmji A das 



where in the MS scheme 



P(a s , e) = f, aS 2 = -ea s + (3(a s ) ■ (6) 
dmjjL z 



Expanding Eq. @ in e and keeping only the leading order term one gets 

7iv(a 5 (/i 2 ) = -— - — r^ } (a 5 ) (7) 
a in as 

where is the coefficient of 1/e in the expansion of T N (as(^ 2 ), e) 

T N (as^),e) = l + Y,^^ 1 - (8) 

e 

In the parton picture, the physical hard scattering cross section ajj with incoming hadron 
H is obtained by convoluting the partonic cross section <j®T with "bare" parton densities 
q B,H . The finite quantity (r ® q B ' H )(x) = qf/ff(x, /i 2 ) is then interpreted as the physical 
(renormalized) number density in the longitudinal momentum fraction x of parton type 
/ in hadron H. This factorization scheme is shown diagrammatically in Figure 1. 
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Figure 1: Factorization scheme 
The parton evolution equation in moment space becomes simply 

rfln 2 = 77v(a s )g7v(^ ) , (9) 

and therefore the anomalous dimensions 7 at (as) defined through Eq. (Q) can be interpreted 
as the moments of the Altarelli-Parisi splitting function 

lN,q/q(a s ) = I dxx N ~ 1 P g/q (x) 



(10) 



where we consider the density of a quark within a quark. These relations allow to calculate 
P q / q (x, as) in a power series in as 

2 



which results in the following expansion of T q / q (x, as, e) 



rg/ g (x,a 5 ,e) = 5(1 -ac) 



«5 

2vr 



+ ... +0 



:i2i 



Here P q / q ( x ) is the first order Altarelli-Parisi splitting function: 

The CFP method for calculating T q / q (x, as, e) in terms of a well-defined projection of 
the 2PI kernel KQ b ' aa ,(k, k') (upper lines are non-amputed, lower lines are amputed) of the 
ladder diagram 



1 



K = K (1-(1-V)K 



(13) 



X-VK 1 

consists in using a projector V = V e <8> V n which can be considered as a pinching of 
the rungs of the ladder. This symbolic representation of the projector V indicates a 
projection acting on the e— dependence (V e ) and a projection acting on the spin, colour 
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and momentum dependence P n of the kernel. The rungs above the pinch are put on-shell 
and considered as external on-shell legs for the upper part of the diagram with averaging 
over spin and colour labels. The rungs below the pinch remain non-amputed. They have 
to be summed over the spin and colour indices in conjunction with the insertion of the 
operator for quarks and — g^ v for gluons and have to be integrated over the off-shell 
loop momentum k. Finally, the projector extracts the single pole term of the whole 
expression, setting e to zero in its residuum according to Eq. (g). The singular non-singlet 
quark density is therefore given by the expression 



r q / q {x,a s , -) = zJ<5(l-x) 



f d m k kn^ 
+ PP i x / TT^ 6 X Trace 



fx K 



} (14) 



Akn 1 - VK 

where PP means to take the pole part. In fixed order perturbation theory we can expand 
r - 1 = K/(l - VK) in terms of the 2PI kernel K as 

r - 1 = VK + V(K ) 2 - V(K VK ) + . . . (15) 

These equations are the basic ones to calculate the anomalous dimensions in higher orders. 
In two-loop order the contributions can be classified into virtual and real contributions. 
The virtual contributions are diagrams where only one internal parton line is cut and 
the diagrams have virtual subdiagrams. The real contributions are defined in terms of 
diagrams where two parton lines are cut and without any virtual subdiagram. In the CFP 
method it is convenient to treat separately the contribution from Zp, the renormalization 
factor of the external leg, and hence to introduce the notation 

T q/q (x,a s ,^) = Z F f q/q (x,a s ,^) 

Fq/q(x,ais,-) = 5(1 -x) 

-\ { (^-\ P®(x) + ! (^-Y P^(x) + 0(al)\ + 0(± 



) P «/> } + 2{Y) p ^ x) + ° {as) + (16) 



z F = + o{ 1 



Pq/ q (x,a s ) = P q / q {x,a s )+{,{a>s)5(l-x). (17) 

An important simplification in the the CFP method is that one has to calculate only 
P^l(x) because the contribution of the wave function renormalization factor Zp at x = 1 
can easily be obtained from fermion number conservation and momentum conservation 
sum rules. 
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3 Principal value versus Mandelstam-Leibbrandt 
prescription 



As we already noted, a crucial ingredient of the CFP method is the use of the light-cone 
gauge (or light-like axial gauge), defined by 

nM*(x) = ; n 2 = (18) 

which formally leads to a gluon propagator of the form 



q z + it qn 



The 1/qn— factor in the gluon propagator gives rise to the so-called "spurious poles", 
singular terms in both the real and virtual contributions. Although these singular con- 
tributions must cancel in gauge invariant quantities, one has to use some regularization 
method for the evaluation of the individual contributions. As yet, in all calculations based 
on the CFP method, the principal value (PV) prescription has been used which is defined 
as 

1 -> lim 1 ( 1 + 1 \- qn (20) 

qn 2 \qn + iS(pn) qn — iS(pn) J (qn) 2 + 5 2 (pn) 2 

where p 2 = ; pn ^ ; pr = . 

It is well-known]] however that the PV prescription is incompatible with Wick rotation 
and hence power counting breaks down. Nevertheless, the spurious poles appearing as 
In 5, In 2 S— terms in individual contributions cancel in the sum of all virtual and all real 
contributions. Another difficulty of the PV method is related to the ultraviolet renormal- 
ization constants which become dependent on longitudinal momentum fractions and on 
In 5. CFP performed the ultraviolet renormalization by subtracting all ultraviolet poles, 
normal and spurious ones, calling this procedure "a phenomenological rule" . 

After the work of CFP it was pointed out that the principal value prescription is 
not consistent with canonical quantization in light-like axial gauge |2l], ^T|. Correctly 
performed canonical quantization leads to the ML prescription defined as 

1 , qn* , , 
► lim 21 

qn 77^0+ qn qn* + in 



3 The subject of quantization and renormalization in noncovariant gauges and their use in perturbation 
theory is described in refs. [|l6|, |l7], |l8| |. 
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where the vector n* is "dual" to n M : 



»fy = ("o, n);n* ll = ( n o, -n) ; (n*) = 0. (22) 

The ML prescription also appears in the literature as 

— -> lim — *; — - -. (23) 

qn v^o + qn + in sign(gn*) 

In the present calculation, these alternative definitions lead to the same result. A crucial 
property of the ML prescription is that the spurious poles are placed in the complex 
go - plane in the same way as the "usual" covariant poles. Therefore, Wick rotation and 
power counting theorems for UV divergences remain valid. The propagator in the ML 
prescription can be decomposed into a term corresponding to the propagation of the 
physical polarizations and into a term which describes the propagation of scalar and 
longitudinal gluons in the qn = plane 

< 0\T{A;(x)A b u (0)}\0 >=< 0|T{T;(:r)T>)}|0 > + < 0|T{L-(a;)^(0)}|0 > (24) 

where 

i5 ab fd A qe iqx 



<0|T{T»T^(0)}|0> = 



and 



(2tt) 4 7 q 2 + ie 

( . (j^nQv + q^) 2qn* n^g Qv + n v g^ 2 \ 

-9V» + — - — 2 1 25 

\ q± nn* qj_ J 



<0|T{ W Lt(0)}|0> = -—j-^- 



( {n^q v + q^n u ) 2qn* n^g Qu + n v g Qll 2 2 \ 
\ q\ nn* q\ Q Q±) J ■ 

Adding up these contributions and using q 2 + q± = 2q + q~ = 2(qn*)(qn) /nn* we get the 
axial-gauge propagator with ML regularization 

KM = 7^ / f-fe + ( "^ + ^"' ) • (27) 

^ (27r) 4 J q 2 + ie\ qnqn* + ie J 

The discontinuity of this propagator can be decomposed into the physical axial-gauge 



contribution and an unphysical contribution 

disc 



= 27rfl( Q b)|-^ + — g 1%) 



<7 2 



|^(^n^)| 2 , 
[ nn* g| J 
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The term proportional to the delta function S(q 2 + q\) is the so-called axial ghost contribu- 
tion. It has been shown by Bassetto et al. that all vectors of the physical Hilbert-space 
are annihilated by the creation operator of these degrees of freedom. Therefore, similarly 
to the Gupta-Bleuler ghosts of QED, they decouple from the S-matrix. The ghosts have 
negative mass squared, indefinite metric, they live in the qn = plane and their polariza- 
tion sum is 

Z<(^M = -K {JM ^ A - (») 
A t^ 2 nn* qi 

In the ML prescription there is a cancellation between the standard axial gauge con- 
tribution and the axial ghost contribution: In the limit q\ — > the spurious poles l/q± 
appearing in the denominator of the polarization sums cancel each other as can be seen 
from Eq. (p8|). As a consequence, when contributions of diagrams with two cut gluon lines 
are calculated, the result must be independent from the regularization scheme of these 
spurious poles. One may use for example dimensional regularization or a PV regulariza- 
tion. The PV regularization can be chosen such that the standard axial gauge contribution 
from the first term in Eq. (|28|) becomes identical to the one of the PV calculation of CFP 
where no ghosts are included at all. Choosing the same PV regularization for the ghost 
part, one finds that the spurious poles cancel within the sum of all ghost diagrams, but a 
finite contribution is remaining from the ghost diagrams. 

Using only dimensional regularization, mixed products of spurious, soft and collinear 
poles are obtained and one cannot see the details of the cancellation mechanism. In 
addition, one deals with poles of higher powers such that the integrals have to be calculated 
in higher order in the regularization parameter e. If we choose PV regularization, however, 
the various (spurious, soft, collinear) singular contributions can be clearly separated from 
each other. In particular, the result we get with ML prescription for the real contributions 
differs from the CFP result only because of the axial ghost contributions. 

One can find a similar decomposition also for the virtual contributions with the help 
of a formula found in refs. KI3, I2BJ . Let us consider the integral 



»...>,-,) = /f V(j+ti) , ' (i+u> (30) 

defined in Appendix A. Using exponential parametrization for the propagator denominator 
factors and the ML prescription to regulate the 1/qn— denominator leads to 



J^ ,ML (ki . . . k n -i) = — da . . . da n -i exp{i V] ajk 2 — iR 2 / z — ze} 

i n Jo 

. f d m expjizq 2 } 

J q+ - R+fz + irj sign (q~ -Rr/z) 
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Ri> = ai kt + ... + a n _i^_! 

Z = do + • • • + O n _ i . 

The q— integral 

exp{iz q 2 } 



JML = J d q + dqS m -Vq ± 



q + — R + / z + irj sign(g~ — R / z) 

has been evaluated using the Cauchy theorem and integrating first over dq + , then over 
dq~ . The integral over the transverse momenta is a standard Gaussian integral, and one 
obtains finally 



Jf h = B+r^' fl - eM^R + R~/z - 2 V \R-\}} (32) 
q R + +mz sign R 1 J 



R + + irjz sig 

leading to 



Jn' ML (h ■ ■ ■ k n-i) = in +f-i J da ... da n _ x exp{i ^ ajk 2 - ze} ■ z 1 2 

exp{-i#/z} , _ exp{2lR+R - /z _ 2r]lR - l} } (33) 



R + + ir] z sign i? 
whereas PV regularization leads to 

771 

Jn ,FY (ki ■ ■ ■ k n -x) = zrmzi / cifflo • • • rfa«-i exp{« V ai ■ - ze} ■ z 1_ ? 

2"^ 2 1 jo — 

•exp{-ii^} ( ^ +)2 ^ w)2 04) 

We can see from Eq. (p3|) that in the ML scheme there is a cancellation between the first 
and the second term as R + — > 0. Therefore, the result will have no spurious poles and 
will be independent from the regularization method of the spurious singularity. Similar to 
the case of the real contributions, it is useful to use PV regularization for the individual 
spurious singularities even within the ML scheme. We emphasize that the application of 
a PV regularization here has nothing to do with the prescription required before carrying 
out the loop integral. Using PV regularization before doing the d m q— integral leads to the 
result ([34]). 

Comparing Eqs. (|33|) and (|34]) shows that the PV result can be identified as a subpart 
of the ML result by dropping the second term in Eq. (|33|) and taking the real part of 
the l/(-R + + it] z signR~) term. Hence the difference between the results for the virtual 
parts in the PV and ML prescription is given completely by the contributions from the 
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second term in Eq. fl33|). It follows that the difference between the sum of all PV- and 
all ML-regulated virtual contributions must be equal to the sum of the real axial ghost 
contributions. That is indeed what we will find as a result of explicit calculation. 

4 Details of the calculation 

4.1 Diagrams contributing to terms with colour factor C\ 

Since every colour structure defines a gauge invariant contribution we can simplify our 
exploratory study of the application of the ML prescription by considering only the terms 
that are proportional to G\. Besides the obvious advantage of reducing the number of 
contributing Feynman diagrams, this colour structure has two additional simplifying fea- 
tures: i) the ultraviolet counter term is vanishing, ii) the sums of the real and virtual 
contributions are separately free from soft and collinear singularities. The axial-ghost 
contributions, however, remain important and therefore the C\ structure gives a good 
opportunity to study these contributions in isolation from other complications. The con- 
tributing Feynman diagrams are listed in Fig. 2. The axial ghosts are denoted by dotted 
lines replacing cut gluon lines. 
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4.2 Virtual corrections 
4.2.1 Quark selfenergy 

The quark selfenergy (see topology (e) in Fig. 2) can be split into two terms 

where S F denotes that part of the expression which stems from the Feynman part (the 
part proportional to —g^v) of the gluon propagator, S A denotes the remaining axial part. 
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E F can be written in terms of simple loop integrals as 

m - 

where 



£F = fCr-TT^- (<* J i W + T"J£(*)) 06) 



Q* = iirlT(l + ()(-&)-'. 

We use the notation ejjv f° r e when the singular contribution comes from the ultraviolet 
region. The axial part E is given by 



Z A = ~9 2 Cfj^ (2k 2 A J£(k)+ A J&(k)>f fi+ h A) (37) 

where 



J iW = / p '?PjV = gw^ + W + W)- (38) 

The actual values of the form factors Pq, P%, P 2 , P3 are different for the PV and the ML 
prescription and can be found in Appendix A. In terms of form factors we obtain 

* = C*^«^*Fo (39) 

E A = -C F j£^^-(2 Ak 2 (Pi(k)-P (k)) + P2(k)(A t> H A t> A)) 

= ~C F j£y^ ^ (2 Ak 2 (Pi(k) - P (k)) + 2 P 2 (k)(Apn- Akn+^ A)) (40) 

Then the contribution of topology (e) to the parton density can be written asf] (see 
Eq. (ED) 



1 Note that a factor of two has to be included for diagrams which are not symmetric. 
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where 

T^(x,k 2 ,e) = C F ^y 2 Tr[^- fzL frf fif ftd^-h) 
= -C 2 F a 2 s (A7rYT(l + e)(-k 2 y e 

■- {P q/q (x, e) [T (e - 1) + 4 (Pi (A;) - P (k)) + - P 2 (k)} 

+ 2P 2 (k)^}. (42) 

Details about the phase space integral are given in Appendix B. 

P q / q (x, e) is the leading order splitting function in m = 4 — 2e dimensions 

P„ /9 (x, e ) = i±^-6(l-x). (43) 

4.2.2 Vertex correction 

The contribution from diagram (d) can also be split into a Feynman part and an axial 
part 

r> = rj )6 + r£ 6 (44) 

which can written in terms of form factors of simple scalar, vector and tensor loop integrals 
given in Appendix A. Inserting this expression into Eq. ([14]) and using a definition for 
similar to the one we used for diagram (e) we obtain for the Feynman part 

T^\x,k 2 ,e) = -C F {C F -^) a 2 s {A^YT{l + e){-k 2 r 



4 - R + (2 - e)R 2 - (1 - e){R 4 + R 5 )} [x + e (1 - a;)] 



x 



+2R 6 (l-e) 2 P q/q (x,e)} (45) 

and for the axial part 

T^\x,k 2 ,e) = -C F (C F -^)a 2 s (4n)T(l + e)(-k 2 )-* 

1 i (P (k)-P 1 (k))P q/q (x,e)-P 2 (k) 1 r 



xv 1 — x 



+(p 2 /k 2 )-* i (P (p) - P!(p) - P a (p)) P ?/9 (x, e) 
" (5 , 1 + 5 , 4 + (2-a;)5 r 2 -25 , o + i2o)P w (a:,e) 



+ - 

x 



-2{R l + R 2 -R G )\±^}\. (46) 
1 — x j J 

The form factors Ri, Pj, Si; i = 0, . . . , 6; j = 0, ... ,2; I = 0, . . . 4 are defined and given in 
Appendix A both for PV and ML prescription. 
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4.2.3 Ultraviolet renormalization 



We use dimensional regularization to regulate both the ultraviolet and infrared singulari- 
ties. Since we do first the loop integrals we must assume that in m — 4e, e is positive. After 
adding the counter terms for the quark selfenergy and vertex one-loop sub diagrams we 
obtain an ultraviolet finite answer. Then we can analytically continue the result to nega- 
tive values of e and we can go on with the evaluation of the Feynman parameter integrals. 
First we should keep p off-shell and set it on-shell only after ultraviolet renormalization. 
It is known that the sum of the two ultraviolet counter terms for the vertex and quark 
selfenergy diagrams in regular gauges vanishes as a consequence of abelian gauge invari- 
ance. Keeping only the ultraviolet pole part of the form factors T ,R e , Pi (see Appendix 



A) and using Eqs. (|4^ . |45|j4q) one can see that the sum of the ultraviolet contributions is 



zero in the ML case, but is non-vanishing in the PV case. In particular one finds 

4t' ML = C 2 F a s 2 —-[-3 Pqq (x) + 2] 
e uv x 

TijT L = C F (C F -^)a s 2 —-[3 Pqq (x)-2] 

l e uv x 

Tfy PV = C F a s 2 —-p qq (x)[-3 + 4I + 4\nx} 
e uv x 

TijY V = C F (C F -^)a s 2 —-p qq (x)[3-4I Q -2\nx] 

z euv x 



Pqq( X ) 

and so 



1 +x 2 
1 — x 



m(e+d),ML n 
J UV — u 

T<f v + g PV = C 2 a s 2 — - Pqq (x)lnx. (47) 

This is an essential difference between the PV and ML schemes. The leftover UV singu- 
larity ( ff~7| ) in the PV scheme is a remainder from the contributions of the spurious poles to 
the UV renormalization constants. In the ML prescription, there are also additional terms 
which are even not proportional to the Born term p qq (x), stemming from the n*— part of 
the virtual integrals, but those terms completely cancel. Therefore it does not matter in 
the ML case whether they are subtracted or not. In the PV case however, the leftover term 
in (|47|) leads to an additional contribution to the anomalous dimension in the following 
way: Subtracting all UV poles, normal and spurious ones, means that in the external loop 
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integral^] over k 2 , T^ e+d >' pv is replaced by its subtracted value 

rp(e+d),PV __ rp( e +d),PV _ ji(?+d),PV 

Note that J"( e + d )' pv contains a factor of (—k 2 )~ e whereas T^ e y d ^' FV does not. Inserting the 
UV subtracted expression (48) into the k 2 — integral (|4l|) and disregarding the finite parts 
of j^ e + d )' PV j s leading to 







. o. _ * . . _ In 3j . . _ In x 
-k 2 )-*P q/q (x,e)2 Pqq (x)2 



\ ^ ( S )2 PlT-PnW + 2 -^[l-x + Pqq (x) In (1 - x)) + 0(e°) ) . ( 19 ) 



The double pole term has to be dropped since f q / q (x, as, e) is defined through simple poles 
in e (see Eq. (|l6l)). But the single pole term gives a contribution to T q / q (x,as,e) which 
is exactly the difference of the virtual contributions (PV-ML)„ irt obtained by using PV 
respectively ML prescription, as can be seen in Tables [I] and |j. This difference will be 
compensated by the ghost diagrams present only in the ML CclSG, clS will be explained in 
Section Ol 



The complete contributions of the virtual diagrams to f q / q (x, e) are listed in Table [I]. 
Virtual diagrams where the cut line is a ghost line would give a contribution at x — 1 
only, so they do not contribute to f q / q {x, e) (see Eqs. (|16D, (|T 



4.3 Real contributions 

In this section we give some details concerning the calculation of the contributions from 
the diagrams of topologies (b),(c),(h) and (i). Topology (i) represents the subtraction 
terms V(K VK ) in Eq. (|TJ) and consists of two Born diagrams linked by an additional 
projection. 



5 We can consider the last integral over k 2 as either infrared or ultraviolet divergent, depending on 
where we put a cutoff. CFP treat its divergence as a collinear infrared singularity, using an ultraviolet 
cutoff Q 2 . Collins and Soper [jlO| instead introduce an infrared cutoff and consider the last integral as an 
ulraviolet divergent integral. 
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Finite part 








epv 


dpy 


(e + gQ^x 


(e + d) p V 


T)„„( x) 


6 


-6 


7 


-7 


o 





i>„„( x) In 2 X 

Pqq\ J 


o 


o 


-2 


2 


o 


o 


Pqq^x) lnxln(l — x) 


o 


o 


-4 


2 


o 


-2 


T)„„(x) Xu.X 


-1 


1 


o 


o 


o 


o 


V r,r,( X) 111 (1 X) 


3 


-3 


3 


-3 








p qq (x) Li 2 (l- x) 


4 


-2 





2 


2 


2 


Pqq( X ) K 2 /3 


-2 


2 


-2 


2 








xj (1 — x) 


2 


-2 














X 





1 


-3 


4 


1 


1 


1 


-4 


4 


3 


-3 








x In x 


-1 


1 


4 


-2 





2 


In a; 


1 


-1 


-4 


2 





-2 


In (1 — x) 


-2 


2 














Spurious poles 


Pqq(x) I In (1 - X) 








-4 


4 








Pqq(x) I lllX 








-4 


4 








Pqq(x) I 




















Pqq(x) h 








4 


-4 








xl 








4 


-4 








h 








-4 


4 








Singular part 


Pqq(x) I /e UV 








4 


-4 








Pqq{x) \nx/e uv 








4 


-2 





2 


Pqq(x)/e UV 


-3 


3 


-3 


3 










2 


-2 















Table 1: The results for the virtual diagrams after UV subtraction. The virtual contri- 
bution to f^(x, e) is obtained as a sum of the finite entries of this table multiplied by 
— h. (ff ) 2 - The singular part shows the ultraviolet counter terms. 
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4.3.1 Standard contribution 

The typical integral we have to evaluate for topologies (b),(c) and (h) can be written as 
&) = \ji^ll^ d{m ~ 2) * f d ^ l 2)-M^(x,h,l 2 ,k,e) (50) 

where <i$(/i,Z 2 ) is the two-body-phase space for the two cut lines and the matrix 
element for topology (a). In terms of light-cone variables q^ = (q + , q~ , q) 

A 2 

d&{l u l 2 ) = J dl+dl+dl^dl 2 d {m - 2) hd {m - 2) l 2 5{p + -k + -If - It) 

.S^-Vfi + h + l 2 ) 5{k- + Zf + l 2 ) 5{l\) 5{l 2 2 ) . (51) 

Introducing dimensionless parameters Zj = lj~/p + , j — 1,2 and using the relation 
q~ = (q 2 + (f)/(2q + ) on gets 

d$(h,l 2 ) = ,~\ m / — — dhdl 2 5(l - x - zi - z 2 ) 
{2n) m J z\ z 2 

h.2 h.2 fi 12 

.5(k+r 1+ r 2 )8(- + - + ^ + ^). 

X X Z\ z 2 

— * 

The external loop also requires the integration over dk. The transverse momentum inte- 
grals can be carried out easily if we introduce new momenta such that i) the 5— function 
will be diagonal, ii) the propagator denominators appearing in the matrix element be- 
come as simple as possible. The appropriate choice is different for different topologies. 
For example, for toplogy (b) and (h) the good variables are 

-> -> -> z 2 -> -> -> -> 

h — hi Z 2 = — - {hi + h 2 ) hih 2 = hih 2 cos 9 . 

1 - Zi 

For the phase space integral we obtain the form 



(27r) m T 2 (l - e) v ' Jo \ x J 

■ f 1 duu- e (l-u)- e ■ ^}\7 € \ r d6(sm9)- 2e (52) 
Jo V n ^\ 2 ~ e ) J° 

— * 

where PS is defined as PS = J dk d<&(li, l 2 ). The azimuthal dependence of the denomina- 
tor is either trivial or it can be cast into the form 1/(1 + A 2 + 2A cos^), such that we need 
the integrals 

I A = ty t} / ^(sin^)- 2e = l 



A/7rr(i - e) Jo 
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r(i-e) r M (sing)- 



0Fr(i-e)7o l + A 2 + 2Acos# 



Ic(\ 2 ) 



F(l,l + e;l-e, A 2 
i F(l,l + e;l-e, 



' A 2 



for A 2 < 1 
for A 2 > 1. 



(53) 



In the case of topology (b) for example one has 



bu 



b 



1 — u 

Zl(l — X — Z\) 



X 



such that the integration range for the u— integration will be split at u — according 
to Eq. (§BJ). 

Carrying out these integrals one should keep the regularization parameter of the PV 
scheme appearing in Eq. (p0|) to regularize the spurious poles. As we explained above, 
the contibutions from diagrams (c), (b Dl1 ), (h Dl1 ) and (i Dl1 ) are the same for the PV and 
the ML scheme. There are no ghosts related to diagram (c) since it does not contain cut 
gluon lines. The individual contributions are shown in Tables |2| and |||. The new feature 
of the ML scheme is that we should also add the contributions from the axial ghost terms 
(see Eq. @ and Fig. 2). 

4.3.2 Ghost contributions 

The diagrams considered in this section are given by topologies b Dl2 , h Dl2 , h° 21 , i Dl2 and 
i° 21 . Topology b° 21 leads to the same result as b Dl2 because the diagrams are symmetric 
under exchange of l\ ^l 2 . Diagrams with two cut ghost lines only give a contribution at 
x — 1, which has been omitted since it can be obtained more easily from fermion number 
conservation (see Section |^). 

The phase space for diagrams with one ghost line and one gluon line cut is given by 
(for details see Appendix |B.2| ) 



The results for the ghost diagrams are given in Tables ^| and |3|. Table [| shows the results 
for both, the real and the virtual part in the PV and ML schemes, where real means 
standard plus ghost contributions in the ML case. 



ost 




(54) 
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no ghosts 


ghosts 


no ghosts 


ghosts 






(c) 


h Dii 






b D 21 


h D 12 








(b + h-i) u 


(b + h- i) gh 


p q g(x) In 2 (1 - x) 








6 


6 


























Pqq{x) III 2 X 


-1 


-1 





-2 


1/2 


1/2 


1 





2 





1 





p qq (x) In a; In (1 — x) 








-2 


-2 








-2 














-2 


Pqq(x) 111 (1 -X) 


4 





-4 























Pqq(x) IB.X 





-3/2 


























Pqq{x) Li 2 (l- X) 





-2 








1 


1 


-2 

















Pqq( X ) n 2 /3 








-2 


-2 




















x In 2 x 








-1/2 

















-1/2 





In 2 x 








-1/2 

















-1/2 





x\yix 


2 


-7/2 


-4 











2 











-2 


2 


\ax 


2 


-7/2 





-2 








-2 











4 


-2 


{l + x)Li 2 (l -x) 








2 


2 




















l-x 





-7 


1 


-2 














3 





1 





-1 


























x 2 In x In (1 — x)/(l — x) 








-4 


-4 




















Spurious poles 


Pqq(x)I 


4 








4 


-2 


-2 











-4 








Pqq{x) I \UX 








-4 


-4 








2 


2 


4 











Pqq(x)I ln(l -X) 








4 


4 








-2 


-2 





-4 








Pqq{x)h 








4 


4 








-2 


-2 


-4 











I ■ x/(l - x) 








-8 


-8 








4 


4 





8 









Table 2: Results for real diagrams. The contribution to f , (x, e) is obtained as a sum of the entries multiplied by — ^ C\ (|^) 







no g 


hosts 




ghosts 


no ghosts 


ghosts 






(c) 


h D " 




b D 12 


b D 21 


h D ^ 








(6 + /i-i)n 


(b + h- i) gh 










-4 


-4 








2 


2 


2 


2 








p q g(x) kix/e 




















2 





2 











p qq {x) In (1 - x)/e 








-4 


-4 


























(1 + x) lnx/e 








-1 


-1 


























(l-x)/e 








2 


2 



























Table 3: Real diagrams, singular parts 





real sum 


diff rea z 


virtual sum 








{b + h- i) M L 


(6 + h — i)p V 


(ML-PV) rea , 


(e + cf) ML 


(e + gOpv 


(ML-PV)„jri 


sum of ghosts 
















only 


p qq (x) ln 2 x 


1 


1 

















p qq {x) lnxln (1 — x) 


-2 





-2 





-2 


2 


-2 


Pqq(x) Li 2 (l- X) 











2 


2 








x In 2 x 


-1/2 


-1/2 

















In 2 x 


-1/2 


-1/2 

















x In a; 





-2 


2 





2 


-2 


2 


In a; 


2 


4 


-2 





-2 


2 


-2 


1 


3 


3 

















X 


-3 


-3 





1 


1 









Table 4: Comparison of real and virtual parts 



5 Discussion and Conclusions 



The reevaluation of the two-loop anomalous dimension in light-like axial gauge is a good 
testing ground to understand the usefulness and reliability of this gauge. The success of 
the calculation of Curci, Furmanski and Petronzio using the principal value prescription 
for the 1/nq factor is based on some heuristic treatment of the ultraviolet renormalization. 
Since the anomalous dimensions are related to the ultraviolet renormalization of bilocal 
operators, it is important to provide a better field theoretical support for the "phenomeno- 
logica rules" found by CFP. Within the PV scheme one gets spurious UV singular terms 
which are subtracted. The Mandelstam-Leibbrandt prescription eliminates this difficulty 
and is in addition consistent with canonical quantization. Therefore it has to be expected 
that the straightforward use of the Feynman rules as given by the ML scheme naturally 
provides a field theoretically correct treatment, leading to a deeper understanding of the 
heuristic CFP rules. This expectation has been justified by our explicit calculation of the 
two-loop anomalous dimensions proportional to the colour factor C\. 

It is instructive to make a detailed comparison of the evaluations in the PV and ML 
schemes. The necessary informations are summarized in Tables 1, 2, 3 and 4. In evaluating 
one-loop insertions in the ML-scheme, the spurious poles cancel within the loop integrals, 
such that the ML result is free from spurious singularities, whereas the PV integrals can 
be recovered as a subpart of the ML integrals, this subpart being plagued by spurious 
singularities (see Eq. ( |3~4|) and Table |l|). After having carried out the momentum integra- 
tion with ML prescription, one can regularize the spurious poles of Feynman-parameter 
integrals with a PV prescription. If all singularities are regularized with dimensional regu- 
larization the answer for the sum remains the same. The sums of all virtual contributions 
are different in the two schemes and the difference is due to the difference in the ultraviolet 
subtraction. 

The differences in the real contributions can be organized according to Eq. (p8|). The 
first term in this equation is the standard axial gauge contribution. The second term 
defines the axial ghost contribution. The results are summarized in Tables ||, and f|. 
The sum of all real contributions is different in the two schemes, the difference being given 
by the axial ghost contribution which compensates the difference found in the virtual 
corrections. We see again that the individual terms in the ML scheme are more regular. 
If we combine the contributions of the diagrams (b Dl1 , b Dl2 , b D21 ) or (h Dn , h Dl2 , h D21 ) 
or (i Dn , i Dl2 , i D ' 21 ) their sum is separately finite. In addition, the sums of the ghost 
and non-ghost diagrams are also separately finite. The fact that the spurious singularities 
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cancel separately in the sum of the virtual and real contributions is a consequence of gauge 
invariance. 

In summary, the evaluation of the two-loop anomalous dimension in the ML scheme is a 
consistent method. The axial ghost contributions are important to get the correct answer. 
It is remarkable that the phenomenological rule of CFP for subtracting all ultraviolet 
contributions (spurious and non-spurious ones) leads to the same additional terms as 
provided by the axial ghost contributions in the ML scheme. 

It would be interesting to see the differences between the two schemes and the higher 
consistency of the ML scheme also for the remaining colour structures. 
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A Virtual integrals 

We define general n-point integrals, containing no axial denominator 1/qn, by 
J^-^(k 1 ...k n „ 1 ) = Jd m q- ( 



q^ 1 ... q^ 



V (q + &i) 2 ...(q + & n -i) 2 
and n-point integrals containing one axial denominator 1/qn by 



J^-^(k 1 ...k n . 1 ) = Jd m q 



q 2 (q + h) 2 . . . (q + k n -i) 2 qn 



I^(h...k n ) - jd m <i {q + ki y;' { l +kn?qn 

In the calculation of the contribution of the diagrams of type (e) and (d) we need two-point 
integrals and three-point integrals. 



A.l Two-point integrals 

First we introduce a general parametrization of the integrals as follows 

Jf(r) = Q r e T ; T = — + 2 



r 



4{r) = -^Po(r) 



rn 



JU T ) = —{Pi(r)r, + P 2 (r)n; + P 3 (r)n,} 



rn 



Ql = ?7r fr(l + e)(-r 2 )- £ . 

The results for the form factors Pq, P±, P 2 and P3 depend on the regularization of the axial 
denominator 1/qn. 

A. 1.1 Form factors for the ML scheme 

If we use ML prescription we get 

P ML (r) = Li 2 (l) - Li 2 (l - Xr ) (55) 

pML (r) = _Xr^X L (56) 
1 Xr 
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P M I(r) = HL(±_ +2 + ^ (57) 
nn* \e uv 1 - Xr J 

Ps ML (r) = —(-l + ^- + -[U 2 {l)-U 2 {l- Xr )]\ (58) 
nn* V 1 - Xr Xr / 



where Xr is defined as 



Xr 



2rnrn* 
nn* r 2 



A. 1.2 Form factors for the PV scheme 

In the case of PV prescription on gets 

P PV (r) = A + !^2_ /l + /o i n ( r +) + l \n 2 {r + ) + Li 2 {l) (59) 
cuv e uv 2 

P[ v (r) = — + 2 (60) 
P 2 PV (r) = (61) 

P 3 PV (r) = + J__j 1 + / oln(r + ) + l ln 2 (r+) _ 4 + jLi2(1) l 

2rn \e uv e uv e uv 2 I 



(62) 



There is no n*— dependence in this case and we used the definitions 



rn 
pn 

i 



u 



u z + o z 







h = f 1 du^^ = - 1 -l^6--U 2 (l) + 0(6). 
Jo u ~\~ o 2 4 

A. 2 Two-point integrals for special momenta k and p 

The parameters for the two-point integrals with the special kinematics used in the calcu- 
lation can be read off from relations ( 55]) to (0). In this case 

p = n* (63) 

k 2 

{p-kf = pk = — ; k 2 T = -k 2 {l-x) ; k 2 < (64) 

k + = x ; p + = 1 (65) 
2 kn kn* 

* k = VTT = x 

nn* k z 

2p+p~ 

X P = h— = 1 tor p 7^ 0, p T = 0. 

p 2 
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Note also that 



lim = -1 



Since the integrals J 2 {p) and J 2 ^{p) only contribute to the UV counterterm, p 2 has to be 
off-shell in this case. In the infrared region, where e < and p 2 = 0, the integrals J 2 {p) 
and J 2fM (p) vanish due to the overall factor (— p 2 )~ e . 

A. 2.1 Form factors in the ML scheme 



P ML (A;) = Li 2 {l)-Li 2 {l-x) 
P™ L {k) = ■ > ' lUJ 



1 — x 



[euv 1 - x J 

H"(k) = {-1 + ^ + I 

itf"(p) = Li a (l) 
P^(p) = 1 
P 2 ^(p) = — + 1 

P 3 Mi (p) = £-{-2 + Li 2 (l)} 
A.2.2 Form factors in the PV scheme 

prv(k) = J^ + ^-I 1 + I \nx + lln 2 x + U 2 (l) 

Pf V (k) = — + 2 
^(*) = 

P[ V (k) = ^-\^ + — - — -h + h\nx+ l -\n 2 x-A + L l2 {l) 
2kn [e uv e uv e uv 2 

P PV (P) = ^--h + Li 2 (l) 
Pf V (p) = — + 2 

p 2 pv (p) = o 
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2pn Veuv tuv J 



A. 3 Three-point integrals 

The three-point integrals will only be given for the special kinematics needed in the calcu- 
lation, defined by the relations ( |B3| ) to QS5D, since a general form depending on parameters 
like Xr cannot be so neatly obtained as it was the case for the two-point integrals. Below 
we list the integrals in terms of form factors. 



{k P y v 

Js(k,p) 



Q k e (-k 2 )- 1 {R lI f + R 2 k fl ) 

-Q k e (-k 2 )- 1 (R^p v + R A k»k v + R 5 {kpY v + k 2 R^ v ) 
k"p v + k u p^ 



pn 

_oi { 

pn 



-k ) (S 1 !^ + S 2 k fl + S^n^ + S 4 n* 



Higher tensor integrals have been eliminated by Passarino-Veltman reduction. 
A. 3.1 Form factors for the Feynman part 



R 


Ri 


R2 


Rs 


R4 


-R5 


Re 


1 -K 2 

e 2 6 


e 2 + e + ^ 6 




Ri + l t +3 


-(i+i) 







A. 3. 2 Axial part form factors with ML prescription 



S, 



MX 



11 1 

— + -lnx- 2I* 2 (1) - 2Li 2 (l - a;) - - In 2 a; 



■Si 



ML 



1 1 1 xlnx xlnx 



e 1 — x 1 — x 



Li 2 (l) + 



2x 
1 — x 



Li 2 {l - x) + - 



1 x 



2 1 -x 



In x 



S ML 



1 lnx 



2 1 ln 2 x 

Lloi L — X) 

2V ; 2 1 -x 



S. 



ML 



e 1 — x 1 — x 
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qML 1 , x ^ nx 

o 4 — h 

e 1 — x 

A. 3. 3 Axial part form factors with PV prescription 

5 py = \ - -/„ + -lnar + Zi - / lnx - 2Li 2 (l) - 2Li 2 (l - x) - \ ln 2 x 

1 xlnx a; 
e 2 e 1 — x 1 — x 
1 lnx Ii 2 (l — x) 



S[ v = ^--^ + -^—Li 2 (l-x)-Li 2 (l) 



s pv 



2 



e 1 — x 1 — x 



p V 1 k 2 fl 1 \nx I hix , . xLi 2 (l-x) 1 2 N 

* 3 = -2^(7 + 7T^- /l + T37-^ 2(1) - i-x + 2 lnx y 
sr = 



B Phase space parametrizations 

B.l Phase space integral for virtual contributions 

The phase space integral needed for the virtual diagrams, where only one gluon line is 
cut, is given by 

r d m k 

PS mrt = 2nzj—5(x-z) 6((p - kf) (66) 

where 

(p _ k f = _ ^Ak 2 ^5{{p- k) 2 ) = x 5{k 2 T + (1 - x)k 2 ) 

X X 



Since the integrand has no angular dependence, the angular integral is trivial here, so 
J d m k = K rn ^ J dk 2 ^ d\k T \ \k T \ m - 3 = K m _ 2 Jdk 2< ^- 1 - dk 2 T (4) ^ 



and 

_ m-2 

K ~ 2 ^ 

is the surface of a (m — 2) dimensional hypersphere. Thus in m = 4 — 2e dimensions, the 
phase space for the virtual diagrams is given by 



r d m k 

PS mrt = 27rzj —5(x-z)5((p-k) 2 ) 



j^- e - A K 2 -2e J dk 2 dk 2 (k 2 T y e x5(k 2 + (1 - x)k 2 ) 
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The upper limit of the d(—k 2 ) integral (k 2 < 0) is denoted by a large momentum scale Q 2 
whose actual value is irrelevant since only the pole part of the A; 2 — integration is needed. 



B.2 Phase space integral for axial ghost contributions 

The two-body-phase space where one of the cut lines, say l 2 , is an axial ghost is goverened 
by the condition l 2 = 

4tt 2 



PS 9 - 



host 



j——^ J dl + dl+d^dl 2 d {m -%d im ~% d im - 2) k T 5(p + -k + -If- If) 

.s(™- 2 \k T + h + r 2 ) 5{k- + zr + 1 2 ) 5(i+) 8(ii) e(i° 2 ) 



Now use 5(l\) = 



1 _ f 2 



to eliminate / x 



substitute l 2 by £ 
Note that f > 

e(i° 2 ) 

Hence 

p gghost 



= 2p + l 



dL 



2p+ 



2 —r u,t 2 

because of the theta function 9{l 2 ) : 
9(l+ + l 2 ) = 9(l 2 ) forZ+ = 



2tt 2 f dz\ f -> -> -> -> -> -> 
dz 2 S(z 2 )5(l — x — z\ — z 2 ) I dkr dl\dl 2 5 (kx + h + l 2 ) 



(27T 



1.2 1.2 12 

di5C- + ^ + i + ^) 

X X 



I 2 

10 X X Zi 

For the diagram (h Dl2 ), it is convenient to eliminate l 2 by using the 5— function for the 
transverse momenta. For diagrams (b Dl2 ) and (h D21 ) it is more convenient to eliminate 

— * 

kx, leading to a simpler form of the denominators. But the final form of the phase space 
of course is the same. So for definiteness, the substitutions suitable for topology (h Dl2 ) 
will be given below. 

The angle 9 and the parameter (3 are defined by 



(67) 



Mi 



kxh cos 9 



1 



Using 

J dl\ dkr = 
leads to 



7T 



2-2e 



r 2 (i- e ) 



J dil{i{)- e dk 2 T {k 2 T )- 



rfi 



\ZkT(\ - e) Jo 



d9(sin9) 



-2e 



PS ghost = F e J ^5(1 - x - Zl )J dii{il)- e dk 2 T {k 2 T y 
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/^(- + - + £+-)• SL € \ r^( sing )~ 2£ (68) 
Jo x x zi v/7rr(|-e)io 



Now substitute 



2tt 



4-2e 



(27r) m P(l -e)' 



fc£ = \k 2 \-u ; li = \k 2 \± '-{l-u)-y 

Ob 



then 



2 /Cy XU 

P ~ 72 



1? (l-x)(l-«)y- 
This is leading to 

pS ghost = F t \k 2 \ 2 - 2e x- 1+e (l -x)~ e P duu-'ft-u) 1 -' C dyy~ e 

Jo Jo 

■ [ dZ5{-(l-u)(l-y)+t). j^n e \ [" d9(ain0)-* (69) 

JO X ^TTL ~~ e ) J° 

Finally the integral over 9 can be transformed to an integral from to 1 by substituting 
w — |(1 + cos#), leading to 

pS ghost = F ^ k 2\2-2e x -l+e( 1 _ x ye f 1 duu -e^_ u y-e f 1 dyy -e 

Jo Jo 
B(i-U-<) l 1 *"*" {1 - W)] ~ i ~'- <70) 
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